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OVERVIEW 

Drag  is  generated  when  a  fluid  moves  through 
vegetation.  The  drag  creates  velocity  gradients 
and  eddies  that  cause  momentum  losses. 
These  losses  are  significant  for  a  wide  range  of 
flow  conditions,  and  existing  techniques  for  the 
prediction  of  resistance  do  not  take  these  into 
account,  leading  to  underpredictions  of 
resistance. 


attributed  to  the  effect  of  viscous  shear  at  the 
boundaries  of  the  channel.  This  realization  led 
Prandtl  (1904)  to  formulate  a  logarithmic 
velocity  relation  for  open  channel  flows.  His 
proposed  logarithmic  distribution  has  been 
shown  to  be  a  good  representation  of  the 
actual  velocity  distribution  in  channels  and 
canals  that  do  not  have  obstructions  in  the 
water  column  (Chow  1959). 


Concepts  of  drag  were  employed  to 
formulate  two  new  resistance  relations  for 
cases  when  dense  vegetation  is  present  in  the 
floodway: 


n  =  K„Ry 3 


Cd  Ad 
2g 


The  assumptions  Prandtl  used  to  derive  the 
logarithmic  velocity  distribution  are  not  valid  for 
vegetated  floodways.  In  addition  to  the  viscous 
shear  at  the  boundaries  of  the  channel,  the 
vegetative  elements  generate  form  drag  losses. 
The  drag  induced  by  the  vegetative  elements 
impedes  the  flow  and  causes  a  deficit  in  the 
local  velocity. 


for  unsubmerged  vegetation,  and: 


(RHS  Eqn\T)Jg 


for  fully  submerged  vegetation.  Values  for 
CdAd  as  a  function  of  vegetation  type  and  flow 
condition  are  presented  in  an  accompanying 
technical  note. 


PHYSICAL  PROCESSES 

The  primary  resistance  to  the  flow  of  water  in 
unrestricted  open  channel  flows  can  be 


Fischenich  (1 996)  found  that  the  net  effect  of 
this  drag  on  the  velocity  profile  was  consistent 
for  all  of  the  data  present  in  the  literature, 
regardless  of  the  vegetation  type  or  flow 
medium.  Figure  1  presents  velocity 
distributions  for  three  flow  conditions: 
unobstructed,  submerged  vegetation,  and 
unsubmerged  vegetation.  The  deviation  from 
the  logarithmic  profile  for  the  two  cases  with 
vegetation  impeding  the  flow  is  clear. 

Figure  2  shows  velocity  plots  for  several 
vegetation  types  and  flow  conditions.  The 
heights  and  velocities  are  normalized  for 
comparative  purposes.  The  profiles  display 
similar  characteristics;  retarded  velocities  end  a 
near-zero  velocity  gradient  within  the  lower 
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portion  of  the  vegetation  canopy,  a  zone  near  logarithmic  distribution  above  the  vegetation, 

the  top  of  the  vegetation  with  a  very  high 
velocity  gradient,  and  an  approximately 


Figure  1.  Velocity  distribution  for  submerged  and  unsubmerged  vegetation.  Velocity 
distribution  represents  vegetation  condition  to  the  left 
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Figure  2.  Normalized  velocity  distributions  for  various  types  of  vegetation  and  flow 
conditions.  The  combined  data  represent  27  analyses  of  grasses,  23  analyses  of  shrubs, 
and  16  analyses  of  trees 
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Because  vegetative  drag  can  have  a  profound 
effect  on  the  velocity  and,  thus,  the  water 
surface  elevation,  any  expression  of  the  flow 
conditions  in  a  vegetated  channel  must  include 
drag.  The  general  relation  for  drag  is: 

Fa~CdAV2  (1) 


where 

Fd=drag  force  (ML/T2) 
p=fluid  density  (M/L3) 

Cd=an  empirical,  dimensionless  drag 
coefficient 

A=area  of  the  obstruction  normal  to  the  flow 

(L2) 

V  =approach  velocity  of  the  fluid  (L/T) 

Fenzl  and  Davis  (1964)  found  that  the  relative 
influence  of  the  vegetation  drag  and  the  soil 
shear  for  unsubmerged  vegetation  is  a 
function  of  the  height  above  the  bed  z  relative 
to  the  height  of  the  vegetation  h  and  two 
heights  z-|  and  Z2  between  z  =  0  and  h  that 
are  a  function  of  the  sediment  gradation  and 
the  vegetation  properties  (Figure  3).  They 
determined  that: 

1)  For  a  range  of  small  depths  of  flow  ( z 
<  z-|  <  h ),  soil  roughness  is  the  predominant 

source  of  hydraulic  resistance.  Under  these 
conditions,  resistance  decreases  with 
increasing  depth  of  flow. 

2)  At  some  greater  flow  depth  ( Z2  <  z  <  h ), 

resistance  becomes  essentially  independent 
of  small  changes  in  soil  roughness  and 
increases  with  increasing  depth  of  flow. 

3)  For  intermediate  depths  ( z-|  <  z  <  Z2 ), 

resistance  is  influenced  by  both  soil  and 
vegetative  roughness. 

Fenzl  and  Davis  found  that  the  ratio  of 
resistance  due  to  soil  roughness  to  total 
resistance  decreased  from  a  value  of  0.5  at  a 
depth  of  flow  of  0.1  ft,  to  0.07  at  a  depth  of  0.2 
ft.  Thus,  the  significance  of  the  soil  properties 
within  the  vegetated  portion  of  the  floodway  is 
significant  only  for  very  shallow  flows,  and 
vegetative  drag  is  the  dominant  source  of 


momentum  loss  under  most  flow  conditions. 


Figure  3.  Parametric  description  for 
Fenzl  and  Davis  (1964)  relations. 

RESISTANCE  FOR 

UNSUBMERGED  VEGETATION 

A  relation  for  unsubmerged  vegetation  can  be 
formulated  from  the  principle  of  conservation  of 
linear  momentum.  Following  a  derivation 
similar  to  that  for  the  de  Saint  Venant  Equation, 
the  sum  of  the  external  forces  in  a  control 
volume  (CV)  is  equated  to  the  rate  of  change  of 
linear  momentum: 


£F  =  ma  =  m(-F-)  (2) 

Considering  only  the  x-component  of  the 
linear  momentum,  the  right  side  of  Equation  2 
can  be  expanded  to  form  Equation  3  as 
follows: 


(  pVAJAx + 

The  external  forces  include  gravity  (Fg), 
pressure  (Fp),  drag  (Fd),  and  friction  (Ff),  for 
which  the  x-component  which  can  be  described 
as: 


dpV\A 

dx 


px-(pqsus)bx 


dM  _  d 
dt  dt 


Fg  =  PgAAxSo  (4) 
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-pg 


'&\l 

dx  J 


Ax 


(5) 


Equating  the  slope  term  in  Equation  1 1  to  the 
slope  term  in  Manning's  Equation,  a  relation  for 
Manning's  n  is  established  as: 


Fd  =  -jCdAdF2AAx 
£ 

(6) 

F  f  =  -pgAAxS  f 

(7) 

where 

Fg  =  external  gravity  force  on  the  CV 
So  =  bed  slope 

Fp  =  external  pressure  force  on  the  CV 
Fd  =  external  drag  force  exerted  by  the 
vegetation  on  the  CV 
Ad  =  SA/AXx  =  vegetation  density  per  unit 
channel  length  L'1 

Ff  =  external  friction  force  due  to  shear  on 
the  boundary 

Sf  =  friction  slope  (i.e.  the  slope  of  the 
momentum  grade  line) 


1 1/2 


n  =  KnRy3 


Cd_Ad 

2g 


(12) 


Alternative  derivations  of  Equation  12  from 
energy  or  shear  considerations  are  possible. 
The  form  varies  depending  on  the  assumptions 
made.  Equation  12  requires  an  estimate  of  the 
drag  coefficient  Cd  and  the  corresponding 
vegetation  area.  These  are  discussed  in 
further  detail  later  in  this  chapter.  This 
equation  is  applicable  only  when  the  vegetation 
is  unsubmerged  because  it  does  not  account 
for  the  momentum  flux  downward  into  the 
canopy  that  occurs  with  overtopping  flows 
(Figure  4).  The  following  section  addresses 
submerged  vegetation  cases. 


Collecting  these  terms  and  rearranging,  the 
left-hand  side  of  Equation  2  gives: 


IF  =  ApgAx\ 


So-Sf 


2g  dx 


(8) 


Using  Equations  3  and  8,  assuming  the 
seepage  inflow  and  the  boundary  shear  are 
negligible,  and  rearranging  yields 
Equation  9: 


AdCd  _  So _ 1  dV  1  dV  J_cfy 

2g  V2  gV2  dt  gV  dx  V2  dx 


Figure  4.  Example  of  flow  through 
unsubmerged  vegetation 


which  is  the  unsteady,  gradually  varied  version 
of  the  de  Saint  Venant  Equation  for  linear 
momentum  replacing  the  boundary  shear  term 
with  a  drag  term.  The  corresponding  steady, 
gradually  varied  equation  is: 

A*C4_So_  1  dV_  1_  dy 

2g  V2  gV  dx  V2  dx 


and  the  steady,  uniform  equation  is: 


2g 


=  So 


(11) 


RESISTANCE  FOR  SUBMERGED 
VEGETATION 

Equation  12  does  not  work  well  when  the 
vegetation  is  fully  submerged  because  it  does 
not  take  into  account  complex  flow  conditions 
above  the  flow  canopy  that  can  increase  or 
decrease  overall  resistance  depending  on  the 
depth  of  flow  and  the  Reynolds  number. 
However,  normalized  velocity  profiles  for  this 
case  are  definitive,  and  present  an  alternative 
approach  to  the  derivation  of  a  resistance 
relation. 
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Meteorologists  and  fluid  mechanicists  involved 
with  wind  power  generation,  soil  erosion 
control,  crop  management,  and  other  related 
activities  have  studied  the  behavior  of  winds 
inside  and  directly  above  forest  canopies  and 
crops.  A  modified  logarithmic  law  describing 
turbulent  velocity  profiles  has  been  proposed 
by  most  investigators  for  situations  when 
stratification  has  only  a  minor  influence,  and 
can  be  summarized  by 
u  1  \(z-d)  1  Equation  13,  below: 

—  =  — In  - - - 

U *  k  *-  Zo  -*  For  (z-d)/z0  >  1  (13) 
where 

u  =  local  velocity 

U.  =  (x/p)1/2  is  the  friction  velocity 


d  =  zero-plane  displacement 
k  =  von  Karman's  dimensionless  shear  layer 
constant 

z  =  height  above  the  ground 
z0  =  surface  roughness 

Figure  5  graphically  presents  the  parameters  in 
Equation  13.  The  z0  term  is  generally  taken  to 
represent  a  mean  value  across  an  uneven 
surface  and  is  often  eliminated  from  the 
numerator.  The  displacement  thickness  d  is 
important  for  tall  roughness  elements  such  as 
long  grasses,  agricultural  crops,  forests,  and 
thick  riparian  vegetation.  In  the  absence  of 
roughness  elements  or  when  they  are 
sufficiently  short,  the  displacement  height  is 
zero.  The  parameters  zD  and  d  are  usually 
determined  from  measured  wind  profiles. 


Figure  5.  Velocity  displacement  method 
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It  is  customary  to  assume  the  von  Karman 
constant  k  =  0.4,  to  specify  displacement  height 
as  a  function  of  the  vegetation  characteristics 
and  to  solve  for  friction  velocity  and  surface 
roughness  height  by  fitting  Equation  13  to 
measured  data.  Surface  roughness  estimates 
have  been  estimated  for  flow  data  obtained 
over  different  agricultural  crops  and  forests.  A 
variance  in  the  results  approaching  an  order  of 
magnitude  is  common,  even  for  flow  over  the 
same  surface.  Experimentalists  frequently  fail 
to  obtain  data  above  the  wake  region  of 
individual  roughness  elements  (z  >  l.5h); 
sometimes  the  data  are  taken  during  non¬ 
neutral  conditions;  and  often  upwind 
nonhomogeneities  distort  the  measured 
profiles.  Table  1  summarizes  various 
investigators'  estimates  of  z0  and  d.  In  Table 
1,  s  is  the  plant  spacing,  S  is  silhouette  area 
(the  total  area  of  the  plant  projected  normal  to 
the  flow),  Cf  is  the  shear  stand  drag  coefficient, 
and  other  terms  are  as  previously  defined. 


Equation  13  has  been  found  to  be  a  good 
expression  of  the  fluid  velocity  profile  for  that 
portion  of  the  profile  above  about  2/3  h. 
However,  it  does  not  accurately  represent  the 
profile  within  the  lower  portions  of  the 
vegetation  canopy.  For  flow  within  the 
vegetation  canopy,  different  profiles  have  been 
proposed  by  meteorologists  using  first-  order 
closure  models  that  specify  an  eddy  diffusivity  K 
and  a  drag  coefficient  Cd  for  constant  foliage 
distribution.  However,  none  of  these 
expressions  are  consistent  with  the  observed 
zero  velocity  gradient  within  the  lower  half  of 
the  vegetation  canopy.  Fischenich  (1996) 
determined  that  the  velocity  profile  within  the 
canopy  could  better  be  approximated  by: 

u_J  cosh  (PZ)V 
uh  {  coshf  fi)  J 


Table  1.  Formulae  for  Roughness  Length  Zp  and  Displacement  Thickness  d 


|  Author 

um*u(i$69) 

h(0.5(s/S)j 

Massman  (1987) 

1.07(h-d)e**(-k(Cf/2)**5) 

h*f(CdLAI) 

;  Merohe y  (1993) 

O.lh  <  zo  <  0.13h 

0.67h  <  d  <  0.75h 

Otterman  (1981) 

0.5h(1-e**(-s/S)) 

Paeschke  (1937) 

h/7.35 

Seginer  (1974) 

(lh/k)exp((-4k**3)/(lhCda))**.33 

h-lh/k 

Sellers  (1965) 

2.5  +  In  h  (h  in  cm) 

Standhii((1969) 

0.64h 
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where 


(17) 


fi  = 


2 CdAdJtUh 
Kh 


(15) 


In  Equations  14  and  15, 
u  =  local  velocity 

uh  =  velocity  at  the  top  of  the  undeflected 
vegetation 

h  =  undeflected  vegetation  height 
z  =  height  above  the  ground 
£  =  z/h 

Cd  =  a  stand  drag  coefficient 

Ad  =  vegetation  area  based  on  density 

Kh  =  eddy  diffusion  coefficient 

Fischenich  (1 996)  investigated  the 
relationships  between  the  lumped  drag-area 
term  CdAd  and  the  z0  and  d  terms  in  Equation 
15.  He  found  that  these  parameters  could  be 
approximated  as  functions  of  CdAd  and  that,  by 
replacing  the  Kh  term  in  Equation  15,  an 
expression  for  the  mean  velocity  could  be 
reduced  to  a  function  of  the  vegetation  height  h 
the  flow  depth  y  and  the  lumped  drag-area  term 
CdAd  (Equation  16): 


" 

V 

2.5 

u* 

y 

]l.26h2 


e<nc^V 


JHCtAi) 


(  -  \ 

\0-5 

y 

-u-95 

Jz+j  In 

h 

h 

dz 

0.13e(CiA^4>1 

/ 

V 


J 


Fischenich’s  general  equation  is  rather  difficult 
to  work  with,  but  can  be  reduced  to  a  simpler 
form  with  some  limiting  assumptions. 
Assuming  the  CdAd  value  remains  below  unity 
and  the  flow  depth  is  not  more  than  twice  the 
height  of  the  vegetation,  estimates  of  within  5 
percent  of  the  actual  velocity  can  be  obtained 
by  solving  the  integrals  in  Equation  16  to  yield: 


_F 
U . 


15 

y 


( 0.23/?3  ^ 

\C*A*  J 


-(1.2  +  2  CdAd)* 


_y-1.15/?-0.051n 


LV 


-0.95 


)- 


By  noting  that  Manning's  n  value  can  also  be 
expressed  in  terms  of  V/U,  the  right-hand  side 
(RHS)  of  Equation  17  can  be  used  to  obtain  a 
value  of  Manning's  n  for  a  particular  flow  depth 
y  as  follows: 


(RHS  Eqn  17  )yfg 

Equation  18  is  applicable  only  in  cases  where 
the  depth  of  the  flow  exceeds  the  height  of  the 
vegetation.  Initial  results  of  studies  by 
Fischenich  (1996)  suggest  that  the  practical 
application  of  Equation  1 8  is  restricted  to  cases 
where  the  depth  of  flow  is  at  least  l.ly,  though 
additional  verification  is  needed.  Because  the 
flow  depth  must  be  specified,  the  equation  can 
be  solved  by  iteration  when  making  normal 
depth  computations  or  can  be  used  to  construct 
a  NV  card  for  HEC-2  analyses.  The  vegetation 
height  and  a  value  for  CdAd  must  be  known  in 
addition  to  the  flow  depth.  Measurements  of 
the  vegetation  height  are  relatively 
straightforward.  Specifying  appropriate  values 
for  Cd  and  A  is  the  subject  of  an  accompanying 
technical  note. 


APPLICABILITY  AND 
LIMITATIONS 

Techniques  described  in  this  technical  note  are 
applicable  to  stream  restoration  and 
assessment  projects  that  include  flood 
conveyance  as  an  objective,  and  where 
riparian  vegetation  is  subject  to  flooding. 

The  range  of  applicability  for  the  two  resistance 
equations  has  not  been  fully  explored. 

However,  the  assumptions  made  in  their 
derivation,  early  applications  to  field  data,  and 
observations  from  laboratory  studies,  suggest 
some  preliminary  limitations.  Applicability  is 
limited  to  steady  uniform  flow  for  the 
unsubmerged  case  (or  at  least  quasi-steady, 
quasi-uniform  flow).  Properties  that  likely 
define  the  equation’s  applicability  include: 
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relative  flow  depth,  Reynolds  number,  soil 
particle  size,  and  vegetation  density. 

Work  by  Fenzl  and  Davis  (1964)  suggests  that 
a  minimum  flow  depth  of  approximately  0.1  h 
is  needed  for  the  assumption  of  negligible 
shear  to  be  accurate.  The  actual  minimum 
depth  should  be  a  function  of  the  soil  particle 
size  and  vegetation  density,  but  this  relation 
cannot  be  defined  currently.  At  sufficiently 
great  depths,  drag  should  be  only  a  small 
component  of  resistance  and  flow  can  be 
assessed  using  the  relation  represented  by 
Equation  12.  Again,  vegetation  density  would 
play  a  role  in  defining  this  depth,  but  a  value  of 
20  h  should  be  a  conservative  estimate. 
Equation  19  is  suggested  for  y>20: 

Equation  19  is  based  upon  the  relative 
roughness  of  sands,  gravels,  and  cobbles.  It  is 

lr  1/6 

n- - - —  (19) 

4^5.75  log (%) 
h 

not  a  continuous  function  with  Equation  18,  but 
either  should  yield  reasonable  values  of 
resistance  for  flow  depths  in  the  range  of  20h. 

Vegetation  density  is  a  critical  determinant  in 
the  degree  of  drag  experienced  and  can  help 
define  the  applicability  of  the  equations.  A 
closely  related  issue  is  the  definition  of 
vegetation  area  for  the  drag  relation.  The  drag 
coefficients  needed  for  the  application  of 
Equations  12  and  18  also  depend  upon  the 
definition  of  vegetation  area,  as  addressed  in 
an  accompanying  technical  note. 
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